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Abstract 

In models with large extra dimensions, particle collisions with center-of-mass energy larger than 
the fundamental gravitational scale can generate nonperturbative gravitational objects. Since 
cosmic rays have been observed with energies above 10 8 TeV, gravitational effects in the TeV 
energy range can, in principle, be observed by ultrahigh energy cosmic ray detectors. We consider 
the interaction of ultrahigh energy neutrinos in the atmosphere and compare extensive air showers 
from TeV black hole formation and fragmentation with standard model processes. Departures 
from the standard model predictions arise in the interaction cross sections and in the multiplicity of 
secondary particles. Large theoretical uncertainties in the black hole cross section weaken attempts 
to constrain TeV gravity based solely on differences between predicted and observed event rates. 
The large multiplicity of secondaries in black hole fragmentation enhances the detectability of TeV 
gravity effects. We simulate TeV black hole air showers using PYTHIA and AIRES, and find 
that black hole-induced air showers are quite distinct from standard model air showers. However, 
the limited amount of information registered by realistic detectors together with large air shower 
fluctuations limit in practice the ability to distinguish TeV gravity events from standard model 
events in a shower by shower case. We discuss possible strategies to optimize the detectability 
of black hole events and propose a few unique signatures that may allow future high statistics 
detectors to separate black hole from standard model events. 
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I. INTRODUCTION 



In models with large extra dimensions, the fundamental scale of gravity may be around 
TeV energies jll 3> The presence of extra dimensions affects both sub- and super- 

Planckian physics. Sub-Planckian physics is affected by the presence of Kaluza-Klein modes 
;hat lead to deviations from standard model (SM) predictions in perturbative processes 

i u n □ 

a 121 IS |9(- Searches for these effects in collider experiments have placed bounds on the 
fundamental Planck scale, M* > 1.3 TeV for two extra dimensions and M± > 0.25 TeV 
for six extra dimensions [ipj. Additionally, submillimeter tests of the gravitational inverse- 
square law constrain M* > 1.6 TeV for n = 2 jll| . 

Super-Planckian physics involves nonperturbative effects, the most striking being the 

[l2| and other gravitational objects [lj, (3 [3, III | 



possible formation of black holes (BHs) [12j and other gravitational objects 



in particle collisions with center-of-mass (CM) energy larger than the fundamental Planck 



scale. (For recent reviews, see Refs. |l8, 



19 



20, 



2ll|.) If the fundamental scale is of the order 



of a few TeV, the products of BH decay could be detected in particle col 



26 



27, 



28 



30 



iders 



and in extensive air showers of ultrahigh energies 

Ultrahigh energy cosmic rays provide a natural beam of particles with primary energies 
up to and above 10 s TeV that can in principle probe TeV scale physics. The dominant com- 
ponent of ultrahigh energy cosmic rays (UHECRs) is believed to be protons 3^| generated 
in extra-galactic sources. UHECR protons naturally generate ultrahigh energy neutrinos as 
they traverse intergalactic space through photo-pion production off the cosmic microwave 



background (CMB) 



33 



34j . The threshold energy for pion production off the CMB induces 



a feature in the UHECR spectrum known as the Greisen, Zatsepin, and Kuzmin (GZK) 



feature 



35 



The flux of neutrino secondaries from the pion production depends on 
he assumed extra-galactic proton injection spectrum and generally peaks around 10 6 TeV 



331 1341 1371 . 38]. These secondary neutrinos are often called GZK or cosmogenic neutri- 
nos. Here we study the characteristics of extensive air showers initiated by ultrahigh energy 
neutrinos and compare the production of BHs in TeV gravity theories with SM interactions. 

Ultrahigh energy neutrinos provide a useful means to test TeV gravity. In some TeV grav- 
ity models, the neutrino- nucleon cross section, ct u n, is greatly enhanced, leading to larger 
numbers of neutrino- induced air shower events. In fact, the lack of observed neutrino air 
showers can be used to place a bound on the neutrino-nucleon cross section that has been 
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translated into constraints on M* comparable to collider limits |28l.l39|. However, the physics 
of BH formation and evolution in TeV gravity theories is highly uncertain and model depen- 
dent. As we discuss below, the cross section of the process can only be roughly estimated. 
While some choices of parameters lead to the enhancement of neutrino-nucleon cross sections 
compared to the SM, others choices give cross sections for BH formation orders of magni- 
tude below the SM case. Furthermore, the evaporation process of BHs generates additional 
uncertainties on the fraction of the primary energy that is left to generate a shower. Even 
if a limit on the neutrino cross-section can be derived from the lack of neutrino-induced air 
showers (for example, if the cosmogenic neutrino flux is better constrained), translating a 
bound on a v ^ into a limit on TeV gravity parameters is highly model dependent. Therefore, 
the identification of quantum gravity effects based solely on neutrino event rates is not very 
effective. 

Here we take a different approach by modeling the detailed characteristics of extensive 
air showers initiated by BH evaporation on a shower-to-shower basis, with the expectation 
that the large multiplicity of secondaries will lead to detectable signatures. We first calcu- 
late the fragmentation of BH and the spectrum of secondaries. The secondary particles are 
then developed with PYTHIA and AIRES j^lj into observable extensive air showers. 
We find that BH-induced air showers generally differ from ordinary air showers. Differences 
in shower maxima reach ~ 200 g cm~ 2 between BH and SM events which could be easily 
detected if the first interaction point of the air showers were either observed or fixed by 
the interaction. Unfortunately, the first interaction point of high energy neutrinos in the 
atmosphere is neither fixed by the interaction nor detectable. Unlike protons, the interaction 
length of neutrinos in air is quite large, thus neutrinos interact with almost equal probability 
at any point in the atmosphere. Moreover, the first interaction point is not directly observed 
since fluorescence experiments can only detect the air shower once billions of particles have 
been generated while ground arrays only observe the air shower as it reaches the ground. 
Shower observables such as the muon content and the rise-depth parameter give indirect sig- 
natures that can distinguish BH and SM events in large statistics experiments that combine 
fluorescence detectors and ground array detectors. 

In addition to differences in the overall characteristics of air showers, BH formation 
produces some unique signatures since the fragmentation secondaries span most particles in 
the SM. In particular, heavy BHs may produce several r-leptons. Multiple r's are unique 
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to BH formation and may be differentiated in future UHECR observatories. 

The paper is organized as follows. In the first part of §11, we review BH formation in 
TeV gravity and the physics of the BH-induced atmospheric events. The aim of this part 
is to fix notations and make the paper self-consistent. In the second part of §11 we first 
focus on the cross section uncertainties, which make the identification of atmospheric BH 
formation based solely on event rates ineffective. Next we discuss the phenomenology of 
BH evaporation, which is the backbone of the air shower simulations. In §111, we describe 
the Monte Carlo that we have developed and used to simulate neutrino-induced air showers 
in the atmosphere. The main results of the paper are contained in §IV, where we show 
the outcome of our simulations and discuss the differences between ordinary air showers 
and BH-induced air showers. In §V, we briefly discuss possible detection techniques for BH 
formation based on r production in BH fragmentation. Finally, we conclude in §VI. 



II. BLACK HOLE PRODUCTION IN TEV GRAVITY 



In models with n extra dimensions the fundamental coupling constant of gravity is the 
(n + 4)-dimensional Newton's constant 



The observed four-dimensional Newton's constant Ga 



M 



pi 



(1) 

6.707 x 1(T 33 TeV~ 2 and 



the (n + 4)-dimensional gravitational constant G„+4 are related by 

G4 = G n+ 4V n 1 , 



(2) 



where V n is the volume of the extra dimensions. If V n 3> M~ n , it follows that -C Mpi. 
For the appropriate choices of n and V n , M* can be of the order of TeV energies such that 
gravity and the electroweak scales coincide. These models provide an attractive solution to 
the hierarchy problem of high-energy physics. 

If gravity becomes strong at the electroweak scale, particle collisions with CM energy 
larger than a TeV can create BHs 12 1, branes and other nonperturbative grav- 

itational objects 3, 16, Il2|. BH formation dominates the gravitational channel if the 
extra-dimensional space is symmetric whereas branes form in asymmetric cases In this 
paper we only consider symmetric compactification and BH production since brane decay is 
even less understood than BH evaporation. 



A. Cross sections 



The static and uncharged BH in (n+4)-dimensions is described by the (ri+4)-dimensional 
Schwarzschild solution 



ds 2 = -R{r)dt 2 + i2(r) _1 dr 2 + r 2 dtt 



2 

n+2 ) 



where 



R(r) 



1- \-^> 

r 



(3) 



(4) 



The Schwarzschild radius r s of the BH is related to the mass Mbh by 



8r(^) / Mbh 

(2 + n) V 



(5) 



At energy scales sufficiently above M*, BH formation is a semiclassical process. Thus 
the cross section can be approximated by an absorptive black disk with radius r s . For 
a Schwarzschild BH the cross section is 



+BH\s;n) 



F(s)nri 



F(s) 



(2 + n) 



2 

n+1 



1 

n + 1 



(6) 



where y/s is the CM energy of the collision, s* = M 2 , and F(s) is a form factor. Since 
Mbh ^ Mtj h follows that r s ~ M" 1 and the cross section (JHJ) must be interpreted at the 
parton level. The total cross section for a neutrino-proton event is obtained by summing 
over partons: 



(?vp^BH{x m \ n) = ^2 / dxqi{x, -Q 2 ) Vij^BHixs; n) 



(7) 



where qi(x,—Q 2 ) are the parton distribution functions (PDFs) [42], — Q 2 is the four- 
momentum transfer squared, x is the fraction of nucleon's momentum carried by the parton, 
and yjsx m = M B H,min is the minimal BH mass for which the semiclassical cross section is 
valid (generally M B H,min ~ few M*). 

Equation (J7J) should be interpreted with care as the total cross section value is affected by 
several sources of uncertainty. The first uncertainty comes from the approximate knowledge 
of the PDFs. For instance, the uncertainty in the gluon distribution (the most uncertain dis- 
tribution) is ~ 15% for x < 0.3 and increases rapidly for large x 43j. Furthermore, the PDFs 



are known only for momentum transfer smaller than 10 TeV. In BH events we expect the 



momentum transfer to be of the order of either the mass or the inverse Schwarzschild radius 



44j . Therefore, the momentum transfer can reach hundreds of TeV in UHECR- induced BH 



events. In the calculation of the total cross section we fix the PDFs for momentum transfers 
above 10 TeV to be equal to the 10 TeV value. Although the dependence of the PDFs on 
the momentum transfer seems quite small (at least for momentum transfers smaller than 10 
TeV), the 10 TeV cutoff on the momentum transfer induces an additional uncertainty in the 
integrated cross section Eq. ©. A conservative estimate of the total uncertainty due to the 
PDFs is ~ 20%. 

A second major source of uncertainty in Eq. derives from the physics of BH formation 
at the parton level, which is presently not well understood. The theoretical uncertainties in 
the dynamics of the process at parton level are parametrized by the form factor F(s) and 



have been summarized in Refs. 



18, 



The two main factors that may affect Eq. (jbj) are the 
uncertainty in the fraction of the initial CM energy that goes into the BH and the presence of 
angular momentum. Numerical simulations for head-on collisions in four dimensions suggest 
that the mass of the BH is smaller than the CM energy of the colliding particles, leading to 
a reduction of the total cross section. Rotating BHs have also smaller cross sections than 
non-rotating BHs. A naive estimate of the corrections due to angular momentum suggests a 
reduction of the cross section of about 40%. On the other hand, the non-relativistic limit of 
two-BH scattering indicates that the geometrical cross section can be enhanced by a factor 
~ 250 — 350%, depending on the spacetime dimension. The classical cross section for photon 
capture can also be used to obtain a crude estimate of the cross section of BH formation, 
suggesting an enhancement of the cross section by a factor ranging from 300% (n = 2) to 
87% (n = 7). 

To our knowledge, all the quantitative results of the past literature are obtained from 



Eq. (J7|) by setting F(s) = 1 and neglecting the PDF uncertainties. (See, e.g., 2^,|2jj,|28|,|29j].) 
This is partially motivated by the fact that an exact estimate of the total uncertainty in the 
BH cross sections due to the combined PDF and parton-level uncertainties is unattainable 
at present. However, the arguments listed above suggest that the cross section uncertainties 
range from ~ 40% to ~ 300% which can significantly affect most results. Throughout this 
paper, we take into account these uncertainties when deriving observables. 

The cross section as a function of M±, M BH)min , E u , and the uncertainties described 
above are shown in Fig.[TJ The upper left, upper right, and lower left panels show <j vv ^bh 
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as a function of Mbh,™,™ for three different incoming neutrino energies [E v = 10 6 , 10 7 , and 
10 8 TeV), for different values of M* (disks for M* = 1 TeV, triangles for M* = 2 TeV, 
stars for M* = 5 TeV, and circles for M* = 10 TeV), and for the case of seven dimensions 
(n = 3, red symbols) and ten dimensions (n = 6, black symbols). The symbols give the cross 
section calculated from Eq. (J7J) setting F(s) = 1 and neglecting the PDF uncertainties. The 
lower right panel shows the cross section as a function of energy for M+ = 1 TeV and ten 
dimensions. The red shaded (inner) region shows the uncertainty in the cross section due to 
the unknown Mbh,™,™ which we vary from M± to 10 M*. The green shaded (outer) region 
shows the uncertainty associated with the BH formation at the parton level and the PDF. 
The solid lines in the graphs show the cross section for SM interactions. The uncertainty 
in the SM cross section due to the unknown PDFs at very small values of x is bracketed by 
dashed lines. 

For a given number of extra dimensions, the total cross section increases with the energy 
of the primary neutrino and decreases when the fundamental scale is increased. At fixed 
energy and M*, the cross section decreases with increasing Mbh,™,™- The overall effect 
makes the cross section at fixed energy vary by many orders of magnitude. For instance, 
the ten-dimensional total cross section at E u = 10 8 TeV spans five orders of magnitude, 
ranging from values of tens of pb to millions of pb, where the lower values are obtained for 
large fundamental scales. At fixed energy and M*, the range in cross section values span 
about an order of magnitude unless Mbh,™,™ becomes comparable to the CM energy of the 
event. In this case, the rate of events is dramatically suppressed and the cross section tends 
to zero. The behavior of the total cross section with the energy at fixed M* is steeper for 
higher values of M B H,min- In the range E u = 10 6 — 10 9 TeV, which is of interest to UHECRs, 
the cross section grows approximately like o vp ^bh ~ E® A ~ X - S , where lower exponents are 
obtained for lower M B H,min- For example, in ten dimensions the behavior of the cross section 
is a up ^BH ~ E° U A1 for M BH ,min = M* = 1 TeV, a up ^ BH ~ E™ 7 for M BH , m in = 5M* = 5 
TeV, and a up ^ BH ~ E° u 77 for M BH , m in = 10 M* = 10 TeV. 

The large uncertainties in the values of <j up ^bh make it quite difficult to discriminate 
between different values of the fundamental scale M* with good precision. The range of pos- 
sible u vp ^BH for a given M* overlaps with the range for larger M* because of the theoretical 
uncertainties. Even if <j up ^ B h were to be constrained by experiments, could not be 
determined unless the degeneracy were removed by additional assumptions on MBH,min an d 



S 



by reducing other theoretical uncertainties in o vv ^bh- The dependence of the cross section 
on n is the least dramatic making it also hard to differentiate between different dimensions. 
In addition, as M* becomes larger than ~ 1 TeV, <j up ->bh becomes smaller than the cross 
section for SM interactions and the probability for BH formation decreases accordingly 



B. BH evaporation products 

Once the CM energy of a neutrino collision with a nucleon in the air reaches the BH 
formation threshold, a BH with initial mass equal to a fraction of the total CM energy may 
form. The distribution of the initial BH masses is given by the differential cross section 

d< dMBH = 2 ' ai ^ BH ^ xs ' n ) • ( 8 ) 

Light BHs are favored over heavy BHs. The typical initial BH mass is usually a few times 
M BH ^ min . Therefore, models with larger (smaller) fundamental Planck scale tend to produce 
higher (lower) mass BHs. The integrated probability of BH formation vs the initial BH 
mass is plotted in Fig.|21 for a neutrino energy E v = I0 7 TeV, n — 6, M+ = 5 TeV and 
MBH,min — 1,3, 5,10 M*. The initial BH mass is very sensitive to the value of M^ min . 
For MBH,min — M+, 90% of the formed BHs have initial mass less ~ 20 TeV, whereas for 
M B H,min = 10 M* the 90% threshold is reached at M BH ~ 80 TeV. 

Once formed, the BH decay phase is expected to proceed in three stages: classical, 
semiclassical, and quantum 22j. In the first stage the BH sheds the hair associated with 
gauge charges and angular momentum. In the second stage the BH evaporates semiclassically 
by emission of thermal Hawking radiation with temperature T#. We assume that most of 



45, \m. The 



the energy is radiated into the brane, as only gravitons can "see" the bulk 
Hawking evaporation ends when the mass of the BH approaches ~ M*. At this point the 
semiclassical description breaks down and the BH may either decay completely by emitting 
a few quanta with energy of order of M* or leave a stable remnant with mass ~ M* j^f. The 
details of this last stage depends on the unknown underlying quantum theory. However, the 
semiclassical decay should lead to most of the observable signatures. During the semiclassical 
evaporation, the BH decays in a time (CM frame) 
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r ~ 



n+3 

l BH " 



(9) 
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FIG. 1: The two upper panels and the lower left panel show o vv ^BuiMBH,minlM^k) for E u = 
10 6 , 10 , andlO 8 TeV. The disks, triangles, stars, and circles are for M+ = 1, 2, 5, and 10 TeV, 
respectively. The red (black) symbols are for n = 3 (n = 6). The shaded regions show the 
uncertainties. The lower right panel shows a up ^BH(E u ) for n = 6 and M* = 1 TeV, with the 
MBH,min range in red (inner shaded region) and the uncertainties at the parton level and PDF in 
green (outer shaded region). The solid lines give the SM cross section, with dashed lines showing 
PDF uncertainties. 
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FIG. 2: Integrated probability of BH formation as a function of the initial BH mass for E v = 10 
TeV, n = 6, M* = 5 TeV, and Mnn^min = 1 (solid line), 3 (dotted line), 5 (dashed line), and 10 
(long-dashed line) M*. 

Assuming a Boltzmann statistics and instantaneous BH evaporation, the BH emits an av- 



erage number of quanta 



23 



(N) 



M, 



BH 



2T, 



H 



(10) 



where the Hawking temperature Th is related to the Schwarzschild radius and to the entropy 
of the BH, S BH , by Q 

TH = n ±1= n ± lM B H_ 



Aixr s n + 2 Sbh 

In TableHJ we list the parameters of typical BHs in ten and six dimensions for different 
choices of M± and of the BH mass (Mbh =5, 10, 50, and 100 TeV). The particle emission 



rate for a BH with temperature Th is given by 

dNi c, \ \ . L 



49 



E 2 



(12) 



dEdt 8tt 2 e E ' T - (-l) 2 ^ ' 

where E is the energy, A c is the optical area of the BH Q], and Tj, q, a Si , and A^j are the 
spin, the degrees of freedom, the greybody factors, and the number of quanta of the particle 
of species i. We neglect particle masses which are generally much smaller than the BH mass. 
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TABLE I: Decay time in the CM frame (r, in units of 10 -26 s), initial temperature (Th in TeV), 
initial entropy (S), average number of produced quanta ((N)), and energy per quantum [Ef (N) in 
TeV) of Schwarzschild BHs for various M* (TeV), Mbh (TeV), n = 6 and n = 3 (in parentheses). 





M BH 


r 


T H S 


(N) 


E/(N) 


1 


5 


0.521 


(0.736) 0.553 (0.282) 8 (14) 


5(9) 


1.11 (0.56) 


1 


10 


1.27 


(2.08) 


0.500 (0.237) 17 (34) 


10 (21) 


1.00 (0.47) 


1 


50 


10.1 


(23.3) 


0.398 (0.159) 110 (252) 


63 (158) 


0.80 (0.32) 


1 


100 


24.5 


(65.8) 


0.360 (0.133) 243 (600) 139 (375) 0.72 (0.27) 


2 


5 


0.107 


(0.130) 


1.22 (0.671) 4 (6) 


2(4) 


2.44 (1.34) 


2 


10 


0.261 


(0.368) 


1.11 (0.564) 8 (14) 


5(9) 


2.21 (1.13) 


2 


50 


2.06 


(4.11) 


0.878 (0.377) 50 (106) 


28 (66) 


1.76 (0.76) 


2 


100 


5.03 


(11.6) 


0.795 (0.317) 110 (252) 


63 (158) 


1.59 (0.64) 


5 


5 


0.013 


(0.013) 


3.48 (2.11) 1 (2) 


1 (1) 


6.95 (4.22) 


5 


10 


0.032 


(0.037) 


3.15 (1.77) 3 (5) 


2(3) 


6.30 (3.55) 


5 


50 


0.254 


(0.416) 


2.50 (1.19) 17 (34) 


10 (21) 


5.01 (2.37) 


5 


100 


0.620 


(1.18) 


2.27 (0.997) 39 (80) 


22 (50) 


4.53 (2.00) 


10 


10 


0.007 


(0.007) 


6.95 (4.22) 1 (2) 


1(1) 


13.9 (8.44) 


10 


50 


0.052 


(0.074) 


5.53 (2.82) 8 (14) 


5(9) 


11.1 (5.64) 


10 


100 


0.127 


(0.208) 


5.01 (2.37) 17 (34) 


10 (21) 


10.0 (4.74) 



Integrating Eq. (|12jl over E gives 
dNi _ . AT 3 



dt 



Ci r «</t-5- ^ r ( 3 )C(3) = CiT Si fi 



C(3)T (n + 3)^(n + l) 



57T 



167T 3 



2 

2~ 



(13) 



where /j = 1 (3/4) for bosons (fermions). Since the observed Hawking emission happens 
on the brane, we use the four-dimensional greybody factors of Ref. j^jj. (Greybody factors 
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54| . See also Ref. |55|, |56( for a 



in n + 4 dimensions have been calculated in Refs. |52l . 
discussion on BH recoil effect.) The values of Sj, q, T s , and /j are listed in TablclTTl 



The number ratio of two particle species i and j is |46| : 



Ni 



c j^Sjfj 

12 



(14) 



Using Eq. (JTUJ), Ni can be expressed as 

Ni = (N) Cir t fi . (15) 

Z^j C j L Sjjj 

The number of each particle species formed for the BHs listed in Table |J for M* = 1 TeV is 
given in Tablc lTTTl 

TABLE II: The values of Sj, q, T Si , fi for the SM particles 



species 


Si 


Ci 


r 

1 SI 


fi 


quark 


1/2 


72 


0.6685 


3/4 


charged lepton 


1/2 


12 


0.6685 


3/4 


neutrino 


1/2 


6 


0.6685 


3/4 


Higgs 





1 


1 


1 


photon 


1 


2 


0.2404 


1 


gluon 


1 


24 


0.2404 


1 


W 


1 


6 


0.2404 


1 


Z 


1 


3 


0.2404 


1 


graviton 


2 


2 


0.0275 


1 



For example, a BH of mass Mbh — 50 TeV and M* = 1 TeV according to Table|l] emits 
63 quanta, each with energy of 0.80 TeV. These quanta are translated into SM particles 
some of which decay or hadronize. The final output of the BH evaporation may contain up 
to ~ 2000 particles. 

III. EXTENSIVE AIR SHOWER SIMULATIONS 

Extensive air showers created by ultrahigh energy interactions in the atmosphere can 
be detected with ground arrays and fluorescence telescopes. Ground arrays record the sig- 
nal which is produced by the particles of the shower reaching the ground. Fluorescence 
telescopes observe the fluorescence light produced by the interaction of the atmospheric ni- 
trogen molecules with the electromagnetic component of the developing air shower. The 
fluorescence method pioneered by the Fly's Eye j^, Q| detector is able to reconstruct the 
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TABLE III: Fragmented number of particle species for n = 6 (n = 3), M* = 1 TeV 



M BH (TeV) 


5 


10 


50 


100 


quark 


3 (6) 


7(14) 


42 (104) 


92 (248) 


c. lepton 


0(1) 


1(2) 


7(17) 


15 (41) 


neutrino 


(0) 


1(1) 


3(9) 


8(21) 


Higgs 


(0) 


0(0) 


1(3) 


3(7) 


photon 


(0) 


0(0) 


1(1) 


1(3) 


gluon 


0(1) 


1(2) 


7(17) 


15 (40) 


W 


(0) 


0(1) 


2(4) 


4(10) 


Z 


(0) 


0(0) 


1(2) 


2(5) 


graviton 


(0) 


0(0) 


0(0) 


0(0) 



longitudinal development of the e + e component of the air shower. Fluorescence detectors 
are currently used by the HiRes and Auger j^] experiments and are planned for the 
future EUSO and OWL Q observatories. This technique provides a good estimate 
of the energy of the primary particle that initiate the air shower, since most of the energy 
of the air shower goes into the observable electromagnetic channel. Another advantage of 
the fluorescence technique is the ability to reconstruct the depth at which the cascade con- 
tains the maximum number of e + e~ pairs, i.e., the depth of shower maximum, X m . This 
parameter is sensitive to the type of primary particle, to its energy and to the interaction 
initiating the cascade. The depth of the first interaction point, X , depends on the total 
cross section of the particle considered. Due to the small values of the neutrino-air cross 
sections, ultrahigh energy neutrinos can induce air showers at any depth in the atmosphere 
such that Xq is arbitrary. 

Neutrinos can interact at any depth in the atmosphere with almost equal probability. The 
interaction length of a neutrino with energy E v = 10 9 TeV is \ v cc — 1.1 x 10 7 g cm -2 for the 
charged current (CC) interactions. This is larger than the column depth of the atmosphere in 
the horizontal direction, which is 3.6 x 10 4 g cm -2 . BH forming interactions do not improve 
this situation as the BH formation cross sections cannot be not much greater than the SM 
values (see Fig.[TJ). For example, if M* = 1 TeV and n = 6, \ v bh — 1.7 x 10 5 g cm -2 . Fig. El 
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FIG. 3: The Xq distribution for 100 neutrinos with E v = 10 TeV interacting in a column depth 
of 3.6 x 10 4 g cm" 2 . The SM CC interaction length is X u cc — 1.1 x 10 7 g cm -2 (solid error bars). 
The BH interaction length is X u bh — 1.7 x 10 5 g cm -2 for M* = 1 TeV and n = 6 (dashed error 
bars). 

shows that the X distribution is flat for SM and BH interactions. Thus the X m distribution 
is also flat. As we discuss below, differences between SM and BH interactions are evident 
in X m — Xq. We can directly compare the values of X m by fixing the value of Xq in the 
simulations. 

We developed a Monte Carlo code to study the air showers induced by BH formation in 
neutrino-air collisions and compare the BH-induced air showers to the SM neutrino-induced 
air showers. The code generates observable secondaries from SM neutrino interactions and 
BH evaporation using the PYTHIA (ver. 6.2) package [40]. These secondaries are then 
injected into the AIRES simulator as primaries for the final air shower. In the AIRES code 
the threshold energy for tracking particles in the air showers are 80 keV for gamma rays, 
80 keV for electrons and positrons, 1 MeV for muons, 500 keV for mesons and 150 MeV for 
nucleons. The geomagnetic field is set to the Pierre Auger Observatory. The "thinning" level 
used in this work is 10~ 5 with a weight limitation of 20. (Thinning is a method commonly 
used in simulations of UHECRs to avoid following the huge number of secondary particles 
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by following only a fraction of them with varying weights. See j4l| for further details.) 



A. SM neutrino-induced air showers 

We simulate the air showers induced by CC and neutral current (NC) interactions by the 
following procedure: 

• The differential cross sections are integrated over the fraction of the total momentum 
of the nucleon carried by the parton (x) for all the possible values of the fraction of 
total energy that goes into the hadronic cascade (y). 

• A value of y is sampled from the previous distribution. The mean value of y at the 
energies relevant for UHECRs is 0.2. 

• The energy of the lepton (CC interaction) or neutrino (NC interaction) in the final 
state is given by (1 — y) E u . The CC lepton is injected into AIRES, r leptons cannot be 
simulated by AIRES. Therefore, we calculate the decay length and use the PYTHIA 
generator to obtain the secondaries, which are then injected at the corresponding 
height at which the r decays. Note, however, that the r particles have a mean energy 
of 0.2E U , so most r's reach the ground without decaying and are not converted into 
observable energy. The NC neutrino is not observable and is not injected in AIRES. 

• The hadronic part of the CC and NC interactions are simulated with PYTHIA. The 
secondary particles are then injected in AIRES. 

B. BH-induced air showers 

A similar Monte Carlo code is used to simulate air showers induced by BH formation. 
The BH simulation follows this procedure: 

• The mass of BH is calculated using the probabilities given by Eq. (JHJ). The gamma 
factor of the BH is 7 = E U /M BH . Different cases of M BH>min are considered. 

• The temperature of the BH, and the energy and total number of quanta emitted in the 
evaporation phase are calculated for different choices of n and M*. The fragmented 
number of particles species is computed (as in TablelTTT|) . 
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• The momentum of each quanta in the rest frame of the BH are calculated assuming an 
isotropic distribution. If the quantum generated is a quark or a gluon, the secondaries 
resulting from the parton cascade of this quantum are calculated using PYTHIA. If 
the quantum is a gauge boson, it is decayed using PYTHIA. The momenta of all the 
particles are then boosted to the laboratory system. 

• The secondaries are injected in the AIRES code to simulate the extensive air shower. 
All the secondaries are injected at the assumed first interaction point except for the 
t particles which are dealt with as in the SM air showers. However, the energy of 
the r particles generated by BH evaporation is generally smaller than the energy of 
the r's generated in the SM process. The decay length of BH r's is comparatively 
shorter than in the SM case. Neutrinos, gravitons, and r's that decay after reaching 
the ground are not observable, thus they and are not injected in AIRES. 

• 7r°'s generated by the hadronization of quarks, gluons, and gauge bosons are imme- 
diately decayed by PYTHIA. This is a good approximation since the average pion 
energy is smaller than the critical energy. Therefore, 7r°'s are more likely to decay 
than interact. 

IV. SIMULATION RESULTS 

A. Neutrino-initiated air showers 

We simulated SM-induced air showers for CC and NC interactions as well as air showers 
from BH production. The showers were chosen to have a zenith angle of 70° and a primary 
neutrino energy E u = 10 7 TeV. The first interaction point was fixed to an altitude of 10 km 
corresponding to a slant depth of 780 g cm" 2 . 

SM neutrino air showers are generally dominated by CC interactions because NC inter- 
actions have lower cross section, a^ c = 0.4 a^ c . Moreover, a large fraction of the primary 
neutrino energy of the NC interaction, (1 — y) E u ~ 0.8 E u , is carried out by the scattered 
neutrino and is not observable. Similarly, the CC scattering produces a high energy 
invisible n that does not contribute to the shower energy. The CC v T interaction produces a 
high energy r that generally does not decay before reaching ground level. For E v = 10 7 TeV, 
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the decay length of the scattered r is ~ 500 km. If the r's were to decay before reaching 
the ground, the air shower would appear as a superposition of showers initiated at different 
heights. (We will return to r decay later in §V.) Therefore, as far as CC interactions are 
concerned, the most easily observed primary is the u e . The secondary electron initiates a 
purely electromagnetic air shower that carries ~ 80% of the primary neutrino energy. These 
showers have similar features to electromagnetic air showers: 

• CC v e air showers are /x-poor. The dominant process for [i production in an electro- 
magnetic cascade is photoproduction. The ratio of the pair production and photopro- 
duction cross sections determines the number of /i's in the air shower. This ratio is 
2.8 x 10~ 3 at 10~ 2 GeV and is expected to be ~ 10~ 2 at ~ 10 7 TeV. 

• CC v e air showers develop slower than hadronic air showers. The number of secondaries 
in pair production or bremsstrahlung interactions is smaller than in hadronic inter- 
actions. Additionally, the Landau- Pomeranchuk-Migdal (LPM) effect js3 , 64, 65, 
also contributes in slowing down the shower development once primary energies reach 
above E LPM ~ 10 7 TeV |6j 

We performed a systematic simulation of SM neutrino-induced air showers and checked 
the characteristics discussed above. Here, we only show the more relevant CC-induced air 
showers and compare them to the BH-induced air showers. The left panel in Fig.|3]shows the 
longitudinal development of CC- and BH-induced air showers. We chose the BH parameters 
MBH,min = 2M + = 2 TeV and n = 6. A difference of ~ 200 g cm -2 in X m is evident 
between the SM and BH events. This large difference results from the combination of the 
large multiplicity and hadronic nature of BH-induced air showers and electromagnetic nature 
of the CC-induced air showers. 

BH-induced air showers generally develop faster than typical SM hadronic air showers 
depending on the initial and minimum masses of the BH that give the number of produced 
quanta. For example, if MBH,min = 2 M* = 2 TeV and n = 6, the average BH mass is 
(M BH ) ~ 7 TeV which produces about seven quanta. If all the quanta are quarks, gluons or 
gauge bosons, the number of secondaries produced is ~ 200. This number of secondaries is 
also close to the mean multiplicity for a SM proton-N 14 collision with energy 10 7 TeV in the 
laboratory frame. However, in SM hadronic interactions most of the momentum is carried 
by the leading baryon, the other 199 particles are softer. In the BH case, the momentum is 
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FIG. 4: Number of e + e vs slant depth for neutrino air showers with E v = 10 7 TeV. The SM 
CC air showers are shown in red (dashed lines) and the BH-induced air showers with n = 6 and 
M BH . 

m i n — — 2 TeV are shown in black (solid lines). The left panel has fixed Xq(CC) — 
X (BH). The right panel has shifted X such that X m (CC) ~ X m (BH). 

equally shared by all the quanta such that the shower produces 200 similar secondaries in 
the first interaction which causes a faster shower development. 

If the first interaction point could be observed, the difference between the BH and the SM 
value of X m — X would be clearly distinguished on an event-by-event basis (see Fig.|Uleft 
panel). However, Xq cannot be directly observed due to limited sensitivity of the detectors. 
The right panel in Fig.H] shows the events with the SM curves shifted by 200 g cm -2 and 
renormalized. The difference between the two cases is no longer apparent. 

CC-induced air showers have large fluctuations in X m — Xq (Fig.HJ). This is mainly due 
to fluctuations in the fraction of primary energy carried by the scattered electron. This 
fraction is usually large such that CC-induced air showers behave often like electromagnetic 
air showers. On the odd occasion that a large fraction of the primary energy is carried 
by partons, the air shower is closer to a hadronic air shower. In addition, the LPM effect 
increases the fluctuations in X m — X of electromagnetic air showers, if the energy of the 
scattered electron (yE u ) is larger than Elpm- On the other hand, the number of particles 
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FIG. 5: The left panels show the X m — X$ distribution for SM and BH air showers with n = 6 and 
MBH,min = 2M* = 2 TeV. The right panels show the distribution of the rise-depth parameter, T, 
for the same showers. The upper (lower) panels correspond to E u = 10 7 (10 8 ) TeV. 

at shower maximum, N max , is proportional to the primary particle energy, which is more 
stable in the electromagnetic cascade case (fluctuations in N max are on the ~ 5% level). 

In contrast, BH-induced air showers have small fluctuations in X m — X and large fluc- 
tuations in N max . The fluctuations in X m — X are consistent with the fluctuations of SM 
hadronic-induced air showers. The large fluctuation in N max (~ 20%) is due to the superspo- 
sition of two effects: (i) Each quantum usually carries a large fraction of the primary energy 
and (ii) some of the produced quanta do not contribute to the shower energy (neutrinos, 
gravitons, /i's, and non-decaying r's). 

The left panels of Fig.|S] show the X m — X distribution for neutrino SM and BH air 
showers. For E v = 10 7 (10 8 ) TeV, the average X m — X for BH-induced air showers is 770 
(840) g cm -2 and 970 (1250) g cm -2 for SM-induced air showers. The spread is 62 (72) g 
cm -2 and 75 (140) g cm -2 . The difference between the BH and CC air showers increases 
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FIG. 6: Scatter plot of T vs N max for CC air showers (red void symbols) and BH air showers (black 
filled symbols) with n = 6 and MBH,min = 2M* = 2 TeV. The disks are for E u = 10 7 , the circles 
for E v = 10 8 TeV, and the triangles for E v = 10 7 5 . 

with the energy because the difference between hadronic and electromagnetic air showers also 
increases with energy (below E v = 10 7 TeV BH and SM showers are indistinguishable). Since 
X is not observable, we define an observable "rise-depth" parameter for each individual 
shower, T = X m — Xq.x, where Xo.i is the slant depth where the shower has 10 % of particle 
content of the shower maximum. T is a more realistic parameter as a discriminator between 
BH- and CC-initiated air showers because it is observable. The right panels of Fig. El show 
the T distributions for SM- and BH-induced air showers at energies E v = 10 7 TeV (upper 
panel) and 10 8 TeV (lower panel). The separation between the distributions is evident at 
E v = 10 8 TeV. This trend is better seen in Fig. |B] where T vs N max is plotted for different 
primary energies. To clearly see the difference in the distributions it would be necessary to 
accumulate a large number of neutrino air showers. For the cosmogenic neutrino flux, we 
would expect one 10 8 TeV neutrino for every dozen 10 7 TeV neutrinos [38] . However, the 
cosmogenic neutrino flux is barely detectable by experiments under construction: at most 
a few events between 10 6 and 10 7 TeV are expected to be detected per year 0, 0, [ll|. 
Either there are larger unexpected fluxes of neutrinos or larger detectors will be needed that 
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can accumulate enough statistics to discriminate between BH and SM interactions through 
the T distribution. 

We also simulated the longitudinal development of /x's for each individual shower. Since 
/x's are detected on the ground, we calculated the /i number for different positions of the 
ground detector relative to X m . In Fig.0 we show the number of /i's vs N max for 50 air 
showers at a depth X m + AX, where AX = 168, 336, and 672 g cm -2 N max is essentially 
proportional to the observed energy. CC-induced air showers are /z-poor because of their 
electromagnetic nature, whereas BH-induced air showers are /i-rich like hadronic air showers. 

To summarize, two features should be used to find evidence of BH formation in extensive 
air showers: the rise depth and the \x content of the air showers. The main differences arise 
from the electromagnetic nature of the CC-induced air showers in contrast to the hadronic 
character of the BH air showers. To take advantage of the differences in [i content and the T 
distribution, an experiment should combine both ground and fluorescence observations for 
each individual air shower. The Pierre Auger Observatory is the first such hybrid detector 
consisting of a ground array which sample the particle content at a given depth, looked over 
by four fluorescence detectors which may determine X m and T. For shower core distances 
larger than ~ 1 km and AX larger than ~ 100 g cm -2 most of the signal recorded in 
the ground detectors is dominated by /i's and thus is directly sensitive to the difference 
in n content between BH and CC-induced air showers. The only challenge for the Auger 
observatory to test the BH hypothesis is the low neutrino flux. If the ultrahigh energy 
neutrino flux is at the level of the expected cosmogenic flux, a larger version of the Auger 
hybrid detector would be needed to test these theories. 

B. Shower dependence on BH parameters 

In the previous sections we compared SM- and BH-induced air showers for n = 6 and 
MBH,min — 2M* = 2 TeV. The previous results can be generalized to different choices of 
these parameters. The two quantities that characterize the BH air showers are the cross 
section and the multiplicity of particles. The cross section uncertainties considered in §11 
affect the first interaction point, but do not affect the shower development. The main factors 
that can change the physical characteristics of the air showers are the multiplicity and the 
nature of secondaries originated from the BH evaporation. 
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FIG. 7: [i number distribution at a depth X m + AX vs N max for 50 simulated air showers with 
E u = 10 7 , 10 , 10 8 TeV. The red void (black filled) symbols correspond to SM-induced (BH- 
induced) air showers (n = 6 and MBH,min = 2M* = 2 TeV). N max gives the observed energy of 
the event. 

The multiplicity is controlled by the mean number of quanta, N q , produced in the BH 
evaporation. Most of these quanta are quarks and gluons that hadronize and initiate a 
number of hadronic cascades with average energy (laboratory frame) E q = jMbh /N q . These 
sub-showers reach a maximum at the same depth. Thus the maximum of the shower, which 
is the sum of all sub-showers, is given by the maximum of a hadronic air shower with energy 
E q . The shower maximum has a logarithmic dependence on the energy: 

E„ 



X m = X + A log 



10 



iV 9 TeV 



+ B. (16) 



The simulations give A ~ 60 g cm -2 and B ~ 311 g cm -2 . A is the change in X m — Xo 
per decade of energy and is analogous to the more commonly used elongation rate, which is 
the change in X m per decade of energy. Our results agree well with the experimental results 
and previous simulations that give an elongation rate ~ 60 g cm -2 [7^ . 
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The number of quanta depends on the BH mass and differs for each individual air shower. 
Moreover, the number of quanta varies with M* and n at fixed energy. The shift in the shower 
maximum of two distinct showers initiated by N ql and N q2 quanta is 

X ml -X m2 =A\og x J^ . (17) 

For instance, if the number of quanta increases by a factor of 3 (10), X m decreases by 29 
(60) g cur 2 . 

The left panel of Fig. |H1 compares the longitudinal development of BH-induced air showers 
for n = 3 and n = 6, with Mbh,™,™ — 2M* = 2 TeV. The primary energy is set to E u = 10 7 
TeV, and Xq is the same for both air showers. The number of quanta produced in the BH 
evaporation decreases for larger n. Approximately three times more quanta are produced 
for n = 3 at fixed Mbh- This translates into a shift in X m of ~ 25 g cm -2 . 

The right panel of Fig-IHI shows X m for M* = 1 TeV and 5 TeV with n = 6 and MBH,min = 
2 M*. At fixed Mbh the number of quanta for M* = 1 TeV is six times larger than for M* = 5 
TeV. However, as M B h is usually slightly larger than M B H,min, the = 5 TeV case starts 
with a more massive BH overall and hence produces a larger number of quanta than = 1 
TeV case (see Fig.EJ). These two effects counteract and compensate each other, leading to 
the same number of quanta for both cases and no shift in X m as shown in Fig.|SJ 

Fig.El shows the variation in the longitudinal development for MBH,min — 2M* and 10M*. 
Mbh changes by a factor five, thus increasing the number of quanta by the same factor. 
This is translated into a shift in X m of ~ 50 g cm -2 , in good agreement with Eq. ()17j). 

The variation of X m with M+ and the number of dimensions have no effect on the con- 
clusions obtained in the previous section. The two parameters discussed to discriminate 
between BH- and SM-induced air showers do not depend on the position of the shower max- 
imum. The observable signatures are based on the difference between the electromagnetic 
nature of the CC-induced air showers and the hadronic nature of the BH-induced air show- 
ers. Therefore, deeply penetrating horizontal hadronic-looking air showers will generally 
signal BH formation. 
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FIG. 8: Number of e + e~ vs slant depth for BH-induced air showers. The left panel shows n = 6 
(black solid lines) and n = 3 (red dashed lines) for fixed MBH,min = 2 M* = 2 TeV. The right panel 
shows M* = 1 TeV (black solid lines) and M* = 5 TeV (red dashed lines) for MsH^min = 2 M* and 
n = 6. 
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FIG. 9: Number of e^e vs slant depth for BH-induced air showers with MBH,min — 

2M* (black 

solid lines) and MBH,min = 5 M* (red dashed lines) at fixed n = 6 and = 1 TeV. 



25 



V. OTHER SIGNATURES OF BH FORMATION 



In the previous sections we discussed the different characteristics of neutrino initiated air 
showers in the atmosphere for SM interactions compared to the formation of TeV BHs. Given 
the uncertainties in the BH formation, evaporation processes, and the inherent fluctuations of 
air showers, clear signals of BH formation are difficult to extract and require a large number 
of neutrino events. The problem is analogous to separating proton-induced air showers from 
gamma-ray air showers, but with an additional unknown Xq. If future experiments can 
observe both T and //'s of a large number of neutrino air showers, a separation between SM 
and BH events could be reached. Given the low expected flux of ultrahigh energy neutrinos, 
hybrid observatories larger than Auger would be necessary. 

As an alternative to a large study of neutrino induced air showers, BH fragmentation 
may be observable via a few events that have no significant background. For instance, the 
production of r leptons in BH evaporation have no significant counterpart in the SM air 
showers. The fragmentation of heavy BHs may allow multiple r production with r energies 
two or more orders of magnitude lower than the primary neutrino energy. This kind of 
process is strongly suppressed in SM interactions. 

One effect of the lower energy of r's produced in BH interactions versus the SM case 
is the shorter decay length of the generated r. As a concrete example, if a neutrino with 
energy E u = 10 7 TeV crosses the Andes mountains towards the Auger Observatory this 
neutrino-induced BH can produce one or more r leptons with energies around 10 5 TeV. 
These r's would decay at a distance of about 5 km from the mountains where the Auger 
Observatory is located. A shower from one such r decay from the direction of the Andes 
would be surprising and even more so if two decays from that direction were to occur. If 
the same neutrino had a SM interaction it could create a single r with about 2 x 10 6 TeV. 
This SM produced r will decay after traveling about 100 km, past the Auger Observatory. 

This example illustrates that for a given neutrino flux and flavor content, the number 
of produced r's may help separate rare events that have a BH origin versus a SM origin. 
Earth-skimming events [Zll lZ3( would also show different energies for the generated r's. A 
significant study of these signatures depends on detailed assumptions of the neutrino flux 
and the detector capabilities and will be more fully addressed elsewhere. 
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VI. CONCLUSION 



We considered the possibility of using UHECR observations to detect effects from TeV 
gravity theories. We focused our attention on the formation and fragmentation of BHs at 
TeV CM energies and found that distinguishing BH formation and SM air showers is much 
more challenging than previously expected. 

The first challenge on this type of study is the unknown details of BH formation and 
fragmentation. The BH formation cross section has large uncertainties and varies by orders 
of magnitude with model parameters that include the number of extra dimensions, the 
energy scale of extra dimensions, and the minimum mass of BHs. In principle, contrasting 
the observed neutrino flux with the expected neutrino flux can help constrain the neutrino 
nucleon cross section, but the uncertainties of the BH cross section limit the translation of 
these constraints into constraints on TeV gravity parameters. 

We showed that BH forming interactions generate very different air showers from SM 
interactions, but the inability of realistic detectors to observe the first interaction point 
hides most of the difference between these air showers. We proposed two parameters that 
show the different characteristics of the two types of air showers: the rise-depth and the 
muon content of the air showers. The BH air showers tend to rise faster, given their large 
multiplicity, and have larger muon contents, given their hadronic nature. A BH air shower 
is similar to a hadronic air shower that can occur at a much higher depth in the atmosphere, 
i.e., a very deeply penetrating hadronic air shower. Deeply penetrating SM air showers are 
dominated by CC processes that generate electromagnetic air showers. SM neutrino air 
showers are similar to deeply penetrating photon-showers. The rise-depth and the muon 
content can help distinguish these characteristics of the SM and BH types of air showers, 
but the process requires a large number of events to overcome the inherent fluctuations 
that generally occur from shower-to-shower. Given that present observatories are not large 
enough to study a large number of neutrino events, the kind of distinction we propose will 
not be achieved in the near future. 

In addition to proposing the study of different air shower characteristics, we suggested 
that unique events can arise from BH formation which are suppressed in SM interactions, 
such as the multiple r generation. The rate for these events is low if the ultrahigh energy 
neutrino flux is at the level of the expected cosmogenic neutrinos. However, unusual air 
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showers from the direction of a mountain chain can signal both a larger flux of neutrinos 
and a departure form the SM interactions. 
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